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Abstract 

We study a minimal model to understand the formation of clusters on surfaces in 
the presence of surface defects. We consider reaction diffusion model in which atoms 
undergoes reactions at the defect centers to form clusters. Volume exclusion between 
particles is introduced through a drift term in the reaction diffusion equation that arises 
due the repulsive force field produced by the clustering atoms. The volume exclusion 
terms can be derived from master equation with a concentration dependent hopping 
rate. Perturbative analysis is performed for both cross-exclusion and self-exclusion one 
dimensional system. For two dimension numerical analysis is performed. We have 
found that the clusterization process slows down due to exclusion. As a result the size 
of the clusters reduces. In this model reaction scheme has algebraic nonlinearity and 
plausible mechanism is also given. 



1 Introduction 

Cluster formation at nanoscales induced by surface defects has been studied extensively in 
recent times ( [ITIo^l5T4TT| and the references therein). It has been found that step edges [US], 
dislocations [5HH] and domain boundaries [HHO] play a very crucial role in cluster formation. 
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In our recent paper it is shown that when Ge is deposited on Si surfaces preferential growth 
occurs at surface defects and domain boundaries [TT]. A reaction diffusion model is proposed 
which qualitatively explains the cluster formation [T2]. In this model surface defects and 
domain boundaries are taken as localized reaction centers in the form of point defects and 
ring defects. Furthermore, simple linear form of the reaction is considered. This can be 
justified under certain approximation. In the studied model clusters and adatoms were 
allowed to diffuse normally with intrinsic diffusion coefficients [TT|[r2]. 

Reaction diffusion models in presence of defects has been studied earlier to model a 
number of phenomena such as trapping of exciton in a crystal at a defect, recombination of 
electron and hole and soliton and antisoliton [13]. In these works trapping reactions have 
been studied in which the reactants get absorbed completely or partially at trapping sites 
(i.e. reaction centers). The authors have focused primarily on the statistical properties 
such as long time behavior and self segregation [T1HT5] . Furthermore, these models describe 
a non-interacting system of particle undergoing reactions in a disordered media. Here we 
would like to emphasize that when we say an "interacting" reaction diffusion system, we 
simply mean interaction between particles other than the reactive interactions (or simply 
reactions). 

Our main aim in this article is two fold. In the first place, we plan to study the effect of 
exclusion in the formation of cluster induced by surface defects. The effect of exclusion in 
a multispecies reaction-diffusion system in the presence of disorder has not been studied so 
widely. This exclusion effect arises due to repulsion among diffusing particles. For the same 
kind of particles, this repulsive effect is incorporated in a mean field way in their diffusion 
coefficients, which is an experimentally determined quantity [39]. But, there is also repulsion 
among dissimilar species. So, this must be taken into account at least in a mean field way 
in any reaction-diffusion system. 

Exclusion effect in homogeneous reaction diffusion systems has been studied by a num- 
ber of authors. In lattice models exclusion is incorporated by restricting the occupancy of a 
site strictly to a single particle. Recently, it is shown that in a lattice system with contact 
interactions there could be discrepancies between the lattice and their corresponding con- 
tinuum model. This arises because in the continuum model the diffusion constant becomes 
dependent on the concentration which may take unphysical values for different lattice types 
and the chosen interaction neighborhood [26]. However, the continuum model agrees well 
for mild to moderate contact interaction strength. 

If we look at normal diffusion the diffusivity is independent of the concentration of the dif- 
fusing particles. However, in a multiparticle system in which there is interaction between the 
diffusing particles the diffusivity can depend on the concentration. In such systems anoma- 
lous diffusion might be observed. It was shown that the critical behaviors of non-equilibrium 
absorbing phase transitions under particle conservation are affected when excluded volume 
interaction is incorporated [2H]- Experimental observations have established that all concen- 
tration dependent diffusion process leads to anomalous diffusion [19l[20]. In an interacting 
multiparticle system, concentration dependent diffusion coefficient appears naturally [21]. 
Nonlinear Fokker-Planck equation has been studied in the past that describes the stochastic 
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motion of a particle in a media whose drift and diffusion terms depends on the probability 
density of the particle pD - fSB] . In the model considered here we have two coupled Smolu- 
chosky equations [3H] in which repulsive force on any one type of particles is generated by 
other species. In the developed Smoluchosky equations the repulsive force on a given type 
of particles is assumed to be generated by the concentration gradients of the other species. 
For a single species system, equation studied here is same in form , developed by other 
authors [231I21]- 

On the same note repulsive interaction between particles can also be seen as an exclusion 
effect as the repulsive force originates from an effective field produced by other particles on 
a tagged particle. Note that any two particle cannot occupy the same position at the same 
time. Hence, this effect can be introduced by a repulsive interaction between the particles 
(i.e. hard core repulsion) as it is done here. We would also like to emphasize here that 
this is an alternative way of incorporating exclusion effects in mean field equations. In this 
article we study a reaction diffusion system in the presence of exclusion and disorder. This 
type of approach has been taken to understand chemotaxis in biological problems [3"0|l3~T]. 
Furthermore, exclusion processes on lattice has been studied extensively in the past to model 
problems in physics, chemistry and biology. It is also shown how these reaction diffusion 
equation can be derived using microscopic principles from the master equation [3l)|l3"T]. 

Another important feature here is that the incorporation of nonlinear cluster formation 
reaction scheme. Since there is no proven reaction scheme for the formation of nanoclusters 
on Si surfaces, we use algebraic nonlinearity in the reaction scheme. The relevant chemical 
kinetic equations are derived in Appendix ( iBl) . 

The organization of the paper is as follows. In section [2] we discuss the theoretical model. 
A perturbative analysis of one-dimensional system is also presented. In Section [3] we study 
the effect of exclusion in a simple diffusion process in the presence of a trap site at the origin. 
We show here that self-exclusion gives rise to concentration dependent diffusion coefficient. 
We draw important conclusions about the formation of clusters in the presence of exclusion 
from this simple set up. Numerical results for both one dimensional and the original model 
of two dimensions are discussed in Section HJ Conclusions are presented in Sectional Further 
scope of this work is also discussed there. 

2 Theoretical model 

We consider a reaction diffusion process on a flat surface on which reaction occurs only in 
the neighborhood of reaction centers. Reaction centers are localized regions on the surface 
where we allow the reaction to take place. Away from the reaction centers there is only 
diffusion. At a reaction center we assume 77 adatom combine to form a cluster. The coupled 
reaction diffusion equation is given by 



d x (D s d x S(x, t) + eS(x, t)<9 x P(x, t)) 





3 



dtP(x,t) = <9 x (P p <9 x P(x,£) + eP(x,t)(9 x ,S(x,t)) 

-/^(x)P(x, t) + K f (x)S(pc, t)\ (2) 

We assume that there is no external flux, J(x, t) — and the initial conditions are given by 
S*(x, 0) = 1 and P(x, 0) = 0. These equation are supplemented by appropriate boundary 
conditions. The diffusion and the drift terms in Eq. (p])and Eq. ([2]) can be derived from the 
master equation (see Appendix [SJ. For e = the process is purely diffusive and describes a 
noninteracting system of particles. When e ^ it describes a system in which particles of 
different species interact via volume exclusion. We have modeled this through an additional 
drift term for each species that depends on the gradient of concentration of the other species. 
This can be pictured in the following manner. Consider an adatom in the vicinity of a cluster. 
Due to thermal noise the diffusion term can be clearly written as D s d^S. It is to be noted 
that in Fickian diffusion arising from the nonuniformity of the chemical potential the form 
of the diffusive term remains same, except that self-diffusion coefficient D s is replaced by 
cooperative diffusion coefficient [39]. In addition to this the adatom experiences a repelling 
force (e > 0) — e<9 x P which appears as an additional drift term in Eq. ([I]) and similarly in 
Eq. (|2J) . We note here that the surface defects help reaction to occur forming clusters in 
its neighborhood. On the other hand the cluster repels adatoms, therby preventing them 
to reach the defect site. So, clearly the formation of cluster involves a competition between 
these two counter processes. The exclusion of one species of particle by the other species of 
particles we call here as cross-exclusion. When exclusion of a particle by their own kind is 
involved we will call it self-exclusion. In our model we have not included the self-exclusion- 
terms as it will only make the diffusion coefficient dependent on concentrations. Later in 
section [3] we will consider a trapping reaction at a static defect to examine the effect of 
self-exclusion. 

Our next aim is to analyze the solution of these coupled equations for a perturbative 
exclusion effect with keeping the reaction scheme linear as it is done in our earlier work |12j . 
We further consider a single point defect at the origin. Note that if e = and rj — 1 the 
above equation becomes linear. We can write Kf(x) = kfS(x) and ifft(x) = fc&5(x). The 
above equation for the linear case with a single defect is exactly solvable and we get. 

S(x,t) = l--^=-Pf s (x,t), (3) 

P(x,t) = -^-P p (x,t), (4) 
where k = {kf/ ' yfD~ a + k b / y / P^)/2 and the function H a (x, t) is given by 

ff„(x,t) = erfc(— - exp(|x|A;/P a + fc 2 t)erfc(^= + kt). (5) 

For the nonlinear case i.e. for finite value of e an analytical solution to the above set of 
equations cannot be obtained in a straight forward way. 
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Let us consider the nonlinear case e 7^ and rj = 1 in one dimension with a single defect 
at the origin. We want to see the effect of a small exclusion (0 < e 1) term. We assume 
that the solution to Eq. ([T]) and Eq. ([2]) can be written as [IU] 

S = So + eS 1 + . . . + e n - l S n ^ + 0(e n ), (6) 

P = P + eP 1 + . . . + e n ~ l P n ^ + 0(e n ), (7) 

where So and Po are solutions of Eq. ([1]) and Eq. (T5]) with e = and is given by Eq. ([3]) and 
Eq. (jl]) in one dimension. We can expand the Eq. ([T|) and Eq. fl2]) in a regular perturbation 
series in powers of e, the resulting equations of order n will be the reaction diffusion equation 
with e = and a source (sink) terms centered at the defect sites that are functions of solution 
of order (n — 1) equations. The general n order equation can be written as 

d t S n (x, t) = D s dlS n (x, t) - K f (x)S n (x, t) + K b (x)P n (x, t) + J s>n (x, t), (8a) 

d t P n {x, t) = D p d?P n {x, t) - K b {x)P n (x, t) + K f {x)S n {x, t) + J p>n (x, t), (8b) 

where J St0 = and J p0 — 0, J Si „,(x, t) and J S) „(x, t) are source (sink) functions that are 
written in terms of SVi-fc, P n -k, ^xSn.-fc and <9 x P„_fc for 1 < k < n. Although these equations 
are linear, solving order by order is still very difficult due to the complicated source (sink) 
terms on the right hand side. It is also not our aim to find a perturbative solution of the 
problem at this point. We can gain ample insight by replacing Eq. (TjQ) and Eq. (T5]) by a 
simpler set of equations. 

For small e we can make the following approximation in Eq. ([1]) and Eq. ([2]) 

ed*S~ed K S , ed^P-ed^Po- (9) 

The resulting equations are a set of linear equations with variable coefficients. However, 
theses equations are still far from being solvable. To simplify it further we shall use the 
properties of the functions So(x,t) and Po(x,t). The gradients, d x So(x,t) and d x Po(x, t) 
are odd functions in x and have a finite discontinuity at the reaction center x = 0. So, we 
have d x S (0-,t) = -d x S (0+,t) and d x P (0~,t) = -d x P (0 + ,t). Also we have \d x S (x,t)\ 
and \d x Po(x, t)\ monotonically decreasing for x G (0, — 00) or x G (0, 00) and as x — > ±00, 
d x S (x, t) = d x P (x, t) = 0. Therefore there exist x = x*, y* > such that < \d x S (x*,t)\ < 
\d x S (0 + ,t)\ and < \d x P (y*,t)\ < \d x P (0 + , t)\ for all t > 0. 

We can make further approximations so that the gradients in Eq. Q can be replaced by 
constants which is valid in some time interval [0, £]. Let us choose x*,y* > and t* G [0,£] 
such that it satisfies H s (x*,t*) = H a (0,t*)/2 and H p (y*,t*) = H p (0,t*)/2 (see Fig. [TJ. 
Solution of these equations give the values of x*and y* at which the concentrations has its 
mean value at some time t* . In our calculations we have taken t* = t/2. Equation Eq. 
can now be written as 

ed^S ~ WpSgn(x), e<9 x P ~ — v s sga(x), (10) 

where v s = \ed^Po(y* ,t*)\ and v p = \ed x ,So(x*,t*)\ and sgn(x) = 2(9(x) — 1/2) accounts for 
the discontinuity at the reaction center x = (9(x) is the Heaviside step function) . Since we 
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Figure 1: Schematic diagram: (a) adatom concentration So(x,t*), (b) cluster concentration 
P (x,t*), (c) exclusion term — ed x S (x, t*) and (d) exclusion term — ed x P (x, t*). Thick hor- 
izontal arrows in (c) denotes the direction of the repulsive force on the cluster particles P 
due to the adatoms. In (d) it denotes the direction of repulsive force on the adatoms S due 
to the clusters. v s = \ed x Po(y* ,t*)\ and v p = \ed x So(x* ,t*)\ are drift speed of adatoms and 
cluster arising due to exclusion . Time t* £ [0, t] and x = x* > and y* > are points 
where the concentrations have the mean value. 



are replacing monotonically decreasing functions by constants the approximation Eq. (jlOp 
is valid in the neighborhood of the reaction center for a small time interval [0,t\. Note that 
in this approximation all the particles are moving into or away from the reaction center at 
constant speeds v s and v p but in the actual case this is not true when the gradients are 
monotonically decreasing. However, this overestimation of v s and v p will not alter the basic 
physics of the problem. We obtain the following simplified reaction diffusion equations 



d t S(x,t) = D s d x S(x, t) — v s sgn(x)d x S(x, t) 

-(k f + 2v s )6(x)S(x, t) + k b S(x)P(x, t), 



d t P{x, t) = D p d 2 x P(x, t) + v p sgn(x)d x P(x, t) 

-(k b - 2v p )6(x)P(x,t) + k f S(x)S(x,t). 



(12) 



Here we note that, two very interesting features arise due to the effect of exclusion. First, 
it gives an extra drift term with velocity which is either into or away from the defect site. 
Secondly, it modifies the reaction rates and the reaction terms become different for two 
reacting species breaking constraints of our kinetic scheme (see Eq. (TfJ and Eq. fl2])). Let us 
define kf = kf + 2v s and kb = k^ — 2v p . The Eq. (ITT]) and Eq. (fl2|) after Laplace transform 
can be written in an abstract notation [32] as the following 

|0( S )> = G s (s) \-K f \<j>(s)) + K b ms)) + \J(s))} , (13) 



G p (s) 



-K b ms)) + K f \<t>(s)) 



(14) 
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where (x\<t>(s)) = (f>(x,s) and (x\ip(s)) = ip(x,s) are the Laplace transform of S(x,t) and 
P(x,t) respectively. The flux term appears due to the initial condition S(x,0) = 1 = 
(x\J(s)). The reaction operators are defined by 5(x — x')k a 8{x) = (x\K a \x') and S(x — 
x')k a 5(x) = (x\K a \x'). The Green's functions G s (s) and G p (s) are defined by 

G s (s) =[s- D s dl + v s sgn(x)d x ] _1 , (15) 

G p (s) =[s- D p 8l - v p sgn{x)d x ] _1 . (16) 

The expressions for the Green's functions in Eq. (j!5p and Eq. ( II 6p are given in O Next 
consider diffusion of adatoms on a surface without defects. If we choose the initial concentra- 
tion S(x, 0) = S(x), we know that it will evolve as Gaussian as there is no cluster formation. 
Now suppose that we introduce a force field at the origin such that it gives rise to a constant 
drift velocity v s in the outward direction. The diffusion of adatoms can be described by the 
following equation 

dtS(x, t) = d x [D s d x - v s sgn(x)d x ] S{x, t) (17) 
In the Laplace domain we can write 

,/ \ G s (x\0) 

1 + 2v s G s (0\0) 

Here (j)(x,s) is the Laplace transform of S(x,t). For brevity the Laplace variable s is kept 
implicit in G s (x\x'). Using the expression for the Green's function Eq. (1431 and E]) in Eq. 
(11 8p and taking inverse Laplace transform, the solution of Eq. ( fTTj) can be written as 

Looking at Eq. ( [19]) we note from the first term that adatoms are pushed away from the 
origin to a distance v s t. The second term has a minimum at the origin and reduces the 
concentration by small amount which is of the order v s . 

Now returning to our reaction diffusion problem we should expect that, in the presence 
of exclusion, adatoms will experience an extra repulsive force which is directed outward from 
the center of the surface defect. From Eq. (1T3|) andEq. (I14p we obtain 

^ x s) = Q(x) _ k,G,{x\«)Q{V) _ (hh ~ k,h)G.(x\0)G r (0\0)Q(0) (M) 

t / G,(»|0)Q(0) [ 



where 



A = 1 + kf + kb + kfh - kfkb (22) 

2D s (-p s + ls ) 2D p (-p p + lp ) AD s D p (- Ps + ls )(-p p + lp y 
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p s = v s /(2D s ), p p = -v p /(2Dp), 7 S = yj p\ + s/D s and 7 P = Jp^ + s/Dp. The function 

Q{x) = (x\G s (s)\J) = 1/s. 

The inverse Laplace transformation of <fi(x,s) and ip(x, s) is performed numerically by 
Talbot method [33] which we denote by Si(x,t) and Pi(x,t) respectively. In Fig. [2] (a) and[2] 
(b) we have plotted the concentrations for the case e = denoted by Sq, Pq, actual numerical 
solution 5*, P of Eq. (CD) and Eq. ([2]) obtained by finite difference method and the solution of 
the modified linear equations Si, Pi. The parameters are e = 0.1, D s = 1, D p = 0.25, kf = 
1.0, kt = 0.1 and t = 1.0. In Fig. Etb), we note that the solution Si and Pi and the numerical 
solution S and P at a point close to the defect site have reduced as compared to the bare 
case, i.e. e = concentrations, Sq, Pq. The current due to diffusion and the drift current 
are in opposite directions for both adatom and cluster which effectively reduces the number 
of particle at the origin. However, this is true only for small value of < e 1. For higher 
values of e the diffusion and the drift currents will be comparable and the higher order terms 
in e will also have a significant contribution (see Fig. [4]). 




Figure 2: (a) Concentrations: Numerical solution S,P, solution without exclusion (e = 0) 
So,Po and solutions of the simplified linear equation Si, Pi, (b) Concentrations close to the 
defect at the origin. 



3 Trapping reaction with self-exclusion 

To understand the effect of self-exclusion we consider a single diffusing species with a trap 
at the origin. We assume that particles interact among themselves, so we introduce a self- 
exclusion term proportional to the density gradient as we have done in the previous section. 
When a particle reaches the trap site it gets trapped with a time-independent rate constant 
k. The reaction diffusion equation for the problem is then 

d t u(x, t) = d x (D d x u(x, t) + eu(x, t)d x u(x, t)) — k5(x)u(x, t). (23) 

We use the initial condition u(x,0) = 1 along with appropriate boundary conditions [35J. 
By rearranging terms in Eq. (1231) the term due to exclusion can simply be absorbed in 
a concentration dependent diffusion coefficient, D(u) = Dq(1 + eu) (here e is redefined as 



S 



e = e/D ). We note that with self-exclusion the current due to diffusion and drift are in the 
same direction as this can be can seen from the expression for the current j(u) = —D(u)d x u. 
For e = 0, Eq. ( I23p has an exact solution [IS]. Let us consider the case of a small reaction 
rate k < 1 and e ^ 0. 

With no loss of generality, setting D = 1 we expand Eq. fl23|) in terms of perturbation 
series in k [ID]. 

U = U + KU\ + k 2 u 2 + 0(k 3 ). (24) 



The solution uq = 1 satisfy the zeroth order equation. The equations for u\ and u 2 are given 
by 

d t u x = (1 + e)dl Ul - 5(x), Ul (x, 0) = 0, (25) 

d t u 2 = (1 + e)dlv,2 - 8(x)ux + ed x (uid x ui),u 2 {x, 0) = 0. (26) 

Equation Eq. (|25p describe diffusion with an external flux which for this case is a negative 
point flux at the origin. In Eq. (T2"6"j) the expression ed x (uid x u\) albeit exactly known is quite 
complicated. It is maximum at the origin and monotonically decreases with increasing values 
of | a; | and vanishes at infinity. Inasmuch as we are interested in the behavior of the solution 
at a finite time, we can replace this term by a point flux at the origin without compromising 
the basic physics. This assumption is valid only for small t. We have 



d x (u\d x u x ) 

where j = lim x ^ d x (uid x u{) = (tt + 2)/ (4n(l 
solving Eq. (1251) and Eq. (1261) we have 



3oS(x), (27) 
e) 2 ). Substituting Eq. (jZZj) in Eq. (ESJ) and 



u(x, t) 



k - en' j 



2%/TTi 

,.2 



=erfc 



\x\ 



2 ^ 7r exp ( 4 (l + e )ty v^Ti \2y/(l + e)t 



4(1 



x 



2(1 + 6) 



+ t I erfc 



t I —x" 
\x\\l — exp 



7T \A(l + e)t 



+ 0(k 3 ). 



(28) 



Let u(x,t) = lim^o u(x, t) be the concentration when there is no self-exclusion. Define 
by Am = u(x,t) — u(x,t), the difference in the concentration. From Eq. ([28]) we note 
that the concentration at the origin has increased due to the flux term j and we have 
An ~ {l-l/^/YT~e)K^tJ^+K 2 ej^t/{TT{l + e))-{l-l/{l + e))K 2 t/A. For e = 0.5, k = 0.2 
at t — 1.0 we have An = 0.019 and the corresponding numerical solutions give Au numerica i = 
0.0154 see Fig. (JHJ). We see that in this case the effect is exactly the opposite (compare Fig. 
[2] for the multispecies case). Furthermore, the width of the depletion zone has increased due 
to the increase in the diffusion coefficient by e. 
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0.90 

-4 -2 2 4 
x 

Figure 3: Concentration u(x,t) numerical solution (A) without exclusion e = 0.0 and with 
exclusion e = 0.5 (B) numerical solution, (C) analytical solution Eq. ( 128]) for n = 0.2 and 
time t = 1.0. The difference between (A) and (C) at x = is Au. 

4 Numerical results 

The well known finite difference method [38] is used to to compute the solution of Eq. ([!]) 
and Eq. ([2]) numerically. In one dimension we will first examine the effect of exclusion and 
nonlinearity on the shape of the concentration profile with a reaction center at the origin. 
The parameters used are D s = 1.0, D p = 0.25, kf = 1.0 and kb = 0.1. In Fig. 0] (a) and 
E(b) we have plotted concentration S(x,t) for rj — 1,2 and 3 at e = 0.2 and time t = 1.0. 
The concentration of S is decreases with increase in 77. The variation of concentrations with 
different values of the exclusion parameter e = 0.0, 0.1, 0.2 and 0.3 are shown is Fig. H] (c) 
and H] (d). Here as we increase e, the concentration S(x,t) decreases. It has already been 
discussed in our study of the modified linear equations (see Eq. ( fTTT) and Eq. ( fl2l) . It is 
shown that exclusion effect modifies the reaction rate at the reaction center and consequently 
the concentration decreases. We also note that change in concentration P(x,t) with e is 
negligible for e C 1. This is also due to the fact that change in P(x,t) due to exclusion 
is not first order in e. We have also found that the width of the concentration profile of 
P reduces as the parameter e is increased. In Fig. (|5]) we have calculated the FWHM 
for the concentration P(x,t) for different values of 77. This clearly indicates that exclusion 
or/and nonlinearity suppress the formation of clusters. In Fig. [6] we have plotted the 
concentration profile of P in two dimensions. Here we have used the same set of parameters 
as in the one dimension case. The number of defects is 100 which is uniformly distributed 
in the region —3 < x < 3 and —3 < y < 3. For a given randomly distributed defects 
we have plotted concentration P(x,y,t). Here also we see that as we go from e = to 
e = 0.2 keeping 77 = 1 fixed the concentration decreases. The concentration is also found to 
decrease as we increase the nonlinearity from 77 = 1 to 77 = 3. In Fig. [7] we have calculated 
the mean concentration averaged over the randomness of defect distribution. The decay of 
concentration S monotonically decreases with time and the rate of its decay slows down as 77 
is increased. Similarly for P its mean concentration increases with time and its concentration 
for any given time t increases with increase in 77. We also note that the mean concentration 
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Figure 4: Variation of concentration with exclusion and nonlinearity. (a) and (b) e = 0.2 and 
r) — 1, 2, 3, (c) and (d) 77 = 1, e = 0.0, 0.1, 0.2, 0.3. The parameters are D s = 1.0, D p = 
0.25, fc/ = 1.0 and k b = 0.1 and t = 1.0. 
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Figure 5: Width of concentration P(x,t) for D s = 1.0, D p = 0.25, kf = 1.0 and kb = 0.1 
and t = 1.0, 77 = 1, 2, 3. 
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Figure 6: Formation of clusters in the presence of exclusion, (ai^) e = 0, rj — 1, (bi^) 
e = 0, rj = 3, (ci )2 ) e = 0.2, rj = 1, (d lj2 ) e = 0.2, 77 = 3 for £> s = 1.0, D p = 0.25, fc/ = 1.0 
and kb = 0.1 and t = 1.0. 



for both S and P decreases with increasing e. This clearly suggests that both exclusion and 
nonlinearity suppress the formation of cluster and the effect of both is additive. 

5 Conclusion 

This model is primarily designed to understand the formation of Ge-cluster on Si-surfaces. 
It is a minimal model for the said type of problems. In this work the effect of exclusion on 
the formation of cluster in the presence of surface defects is studied. In our reaction diffusion 
model we have introduced exclusion through a repulsive interaction between the particles 
of dissimilar nature. This repulsive force is taken to be proportional to the gradient of 
concentration and a set of coupled Smoluchowsky equations is obtained. In the perturbative 
regime a set of modified linear reaction-diffusion equation as an approximation to the actual 
process is considered and this gave us important understanding of the effect of exclusion on 
the formation of clusters in the reaction diffusion processes. The solutions to these equation 
are obtained using Green's function method [33]. The most important conclusion of this 
work is that, both exclusion and nonlinear reaction processes considered here suppress the 
formation of cluster. The effect of self-exclusion on diffusion in one dimension with a trapping 
reaction at the origin is studied. It is found that self-exclusion can give rise to a concentration 
dependent diffusion coefficient as obtained in earlier works (22]. The width of the depletion 
zone is found to increase by e (strength of exclusion). Our numerical studies in one dimension 
showed that the width of the cluster concentration profile decreases with increasing e. In 
two dimensions, the mean concentrations averaged over the surface disorder are calculated. 
It is found that the mean concentration of adatoms S decreases monotonically where as the 
mean concentration of P increases monotonically in time respectively. However, for higher 
exclusion and nonlinearity these mean concentrations are found to decrease with the increase 
in exclusion and nonlinearity in the reaction scheme. 
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Figure 7: Mean concentrations of adatoms S for (a)rj = 1, (b)r/ = 2, (c)^ = 3 and cluster 
P for (d)?7 = 1, (e)rj = 2, (f)r/ = 3 and D s = 1.0, D p = 0.25, fc/ = 1.0 and k b = 0.1 and 
0.25 < t < 1.5. 



There can be quite a few extension of the present work. For example, in this model, 
the repulsive potential is considered to be simple linear function of concentration. But, it 
is indeed possible for this potential to depend nonlinearly on concentration. This needs to 
be explored. Furthermore, we have considered here a simple algebraic nonlinearity in the 
reaction scheme. This is due to the lack of sufficient knowledge of the reaction scheme for 
the formation of Ge-clusters on Si surfaces. As this model can be used in other physical sit- 
uations, different types of nonlinear exclusion potential and reaction scheme can be studied. 
These works are in progress. 



A Derivation of exclusion term from master equation 

The diffusion and the exclusion terms in Eq. ([I]) and Eq. from the master equation are 
derived in the following way. Let us discretize the continuous space into cells of size Ax and 
denote by rii > and rrtj > for all i — 1, 2, . . . , N the number of S and P particles in the i-th 
cell. The configuration of the system can be described by the N- vectors n = (ni, . . . , n^) 
and m = (m 1 ,m 2 , . . . ,m N ). Let N be the set of natural numbers, we define the operator 
Hf : -»■ by 

Hfinx, rij_i, Hi, n i+1 , n N ) = (ni, . . . , n^i, - 1, n i+1 + 1, . . . , n N ), 
H'ini, . . . , n»_i, rii, n i+1 , . . . , n N ) = (m, . . . , n^i + 1, n f - 1, n i+h . . . , n N ), 

(29) 

for alH = 1, 2, . . . , N. Note that for % = 1, N we use periodic boundary condition [3T] . 

We define density dependent hopping rates W^(i) and W^(i) of a particle at the ith site 
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for S and P respectively by 

Wf{i) = w s (l - yL s {m i±1 - mi)), 

W±(i) = wpil-iipim^-m)), (30) 

where w s , w p , /i s and \i v are constants and the superscript ± denotes hopping to the site z±l. 
Here we note that a particle S (say) at cell i has lower hopping rate if the cell i ± 1 contains 
more number of particles P than that in the cell i. Similarly for P this is exactly what should 
be happening when there is an effect of exclusion. Let p(n, m, t) be the probability of finding 
the system in the configuration n, m, at time t The master equation is given by [361137] 

N 

d t p(n, m, t) = J2( n i + X ) \W7(j)p{H+_ x n, m, t) + W+(j)p(Hr +1 n, m, t)] 

3=1 

+( mj + 1) [W-(j)p(n, H+^m, t) + W+(j)p(n, H~ +1 m, t)] 
~ [nj [W s -(j) + W+(j)} + m 3 [W~(j) + W p +(j)])p(n,m,t) 

(31) 

Let Si = J2n n iP( n i m ) t) and Pi = J2n m iP( n , m ; t) where represents sum over all 
configurations having m, . . . , tin, mi, . . . , m^, be the mean number of particles in the ith 
cell. Let us assume that there are no correlation between particles so that we can write 
X)n niUjpin, m, t) = SiSj, J2 Q m i m j p(n,m,t) = PiPj and rumjpija., m, t) = SiPj. Mul- 
tiplying rii through Eq. (I3T1) and summing over all configurations we obtain 

d t Si = w s (S i+ i + Si- 1 -2Si)-w s p JS [Si +1 {P i -Pi +1 ) + Si- 1 {P i -P i - 1 ) 

-S i (P i+1 + P i - 1 -2P i )} (32) 

Now we can substitute Si = S(x,t),Pi = P(x,t) and expand Si±i and Pi±\ in Taylor series 

Si±i ~ S{x, t) ± Axd x S(x, t) + Ax 2 dlS(x, t), 

P i±l ~ P(x, t) ± Axd x P{x, t) + Ax 2 d 2 x P(x, t), (33) 

and ignoring 0(Ax 3 ) and higher order terms, we obtain from Eq. fl32|) 

d t S(x, t) = w s Ax 2 d 2 x S(x, t) - w s fi s Ax 2 [- (2d x S(x, t)d x P(x, t) 

+S(x, t)d 2 x P(x, t)) - S(x, t)d 2 x P{x, t)] (34) 

Now rearranging terms in Eq. (|34l) we obtain 

d t S(x, t) = D s d 2 x S(x, t) + e s (d x S(x, t)d x P(x, t) + S(x, t)d 2 x P(x, t)), (35) 

where D s = w s Ax 2 and e s = 2w s fi s Ax 2 . Similarly multiplying through Eq. (!3~i~|) and 
summing over all configurations we obtain the equation for P(x,t). A point to note here is 
that < /i SjP 1, and for simplicity we have taken e s = e p = e. Using expression Eq. (135*1) . 
the general expression in arbitrary dimension can be written as 

d t S(x, t) = d x [D s d x S(x, t) + e s S{x, t)d x P(x, t)] . (36) 
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B Derivation of kinetic equations 

The reaction scheme is described in Fig El We have 

d[S N ] 



dt 

and 

d[S n ] 



k N -i[S][S N -i]-k N [S N ]. (37) 



dt 



fen-l^n-l] - (A;_ (n _i) + k n [S]) [S n ] + fc_ n [5„+i], Tl = 2, . . . , (JV - 1) . (38) 



We us assume that all intermediate states are in equilibrium. This assumption implies that 
intermediate reactions are very fast. With this assumption we get 

k( n -i)[S] [S(n-i)] - *;_(„-!)[£„] = 0,n = 2, . . . , N - 2, (39) 
[SJ = Tr~~~ [^1 [<Vi] = #n[£] [S„-i], n = 2, 3, . . . , N - 2 (40) 

fc-(n-l) 

where if n = k n /k_ n . 

S + S — S 2 ■■■ S N ^! + S — S N 

^-1 ^-(N-l) 

Figure 8: Reaction scheme. 



Using Eq. (J40]) and Eq. flU]) we obtain 

kcius.[Sf -k N [S N }. 



d[SN\ folN 



dt 

In another approximation we replace [Sn] by J2n=2 S n — P to get the equation. 

C Green's functions 

The Green's function in Eq. (1151) and Eq. (1161) has the following form 

G a (s) =[s- D a 8l + u Q sgn(x)d x ] 1 (42) 

exp(— p a (x — x'))A a (x, x') if x < x' < 

G a (x\x') = { exp(-p a (x - x'))A a (x', x) if x' < x < 0, (43) 

2D a (- 1 p a +-y a ) eX P(Pa( X + X ') ~ 7a ~ X ')) if X ' < < X > 
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and 

{ 2Da (_p ce +7 ce ) eX P(~P"( X + X ') ~ 1<*( X> ~ X )) if X < < 

exp(p a (x - x'))B a (x, x') if < x < x' (44) 

exp(p a (x — x'))B a (x', x) if < x' < x, 

where A a (x, x') = (exp(j a (x - x')) - ^^exp(7 a (x + x')) \ and B a (x, x') = (exp 

Here p a = v a /(2D a ) and 7 a = y^o| + s/D Q are constants. 
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